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An experimental investigation was conducted to examine acoustic receptivity and sub- 
sequent boundary-layer instability evolution for a Blasius boundary layer formed on a fiat 
plate in the presence of two-dimensional (2-D) and oblique (3-D) surface waviness. The 
effect of the non-localized surface roughness geometry and acoustic wave amplitude on 
the receptivity process was explored. The surface roughness had a well defined wavenum- 
ber spectrum with fundamental wavenumber k w . A planar downstream-traveling acoustic 
wave was created to temporally excite the flow near the resonance frequency of an unstable 
eigenmode corresponding to kt* = k w * The range of acoustic forcing levels, f, and rough- 
ness heights, A/i, examined resulted in a linear dependence of receptivity coefficients; 
however, the larger values of the forcing combination t Ah resulted in subsequent non- 
linear development of the Tollmien-Schlichting (T-S) wave. This study provided the first 
experimental evidence of a marked increase in the receptivity coefficient with increasing 
obliqueness of the surface waviness in excellent agreement with theory. Detuning of the 
2-D and oblique disturbances was investigated by varying the streamwise wall-roughness 
wavenumber a w and measuring the T-S response. For the configuration where laminar- 
to-turbulent breakdown occurred, the breakdown process was found to be dominated by 
energy at the fundamental and harmonic frequencies, indicative of K-type breakdown. 


Nomenclature 

C receptivity coefficient, u l{f j/ii a c- 

c p h phase speed of boundary-layer disturbance. 

F nondimensional frequency, 2 Trfu/U^. 

f frequency in Hz. 

h normalized step height, A h n /Ah\ (n — 1,2,3). 

k wavenumber, (a 2 p 2 ) 1 ^ 2 . 

N amplification factor, f* — a,(x; f)dx. 

R Reynolds number base on S r , yJUryc>x v /v. 

R( i Reynolds number base on length scale /, U^l/v. 
U mean streamwise velocity. 

u rms streamwise fluctuating velocity. 

x streamwise distance from model leading edge. 
x v streamwise distance from virtual origin. 
y wall-normal distance from model surface, 
a streamwise wavenumber, 

o, spatial amplification rate. 

p spanwise wavenumber. 

Ah surface roughness step height. 

S r nondimensional length scale, x v /R. 

5* displacement thickness. 
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e normalized acoustic wave amplitude, u ac /U 
<p phase of disturbance velocity. 

7 coherence function. 

7 ] normalized wall-normal coordinate, (y/. x v ) ■ R. 

A wavelength, 2ir/h. 

ly kinematic viscosity. 

a detuning parameter, (a U7 — 

i/> wave angle. 

Subscripts 

ac freest ream acoust ic field. 

/, II lower and upper branch locations. 
t total disturbance component. 

tr laminar- to- turbulent, transition. 

ts T-S wave component. 

w surface roughness condition, 

oo freest ream condition . 

Introduction 

The transition of a boundary layer from a laminar 
to a turbulent state plays an important role in many 
fluid mechanics problems since it ultimately affects 
basic quantities such as heat transfer and skin fric- 
tion. This transition process is affected by many fac- 
tors such as freest ream (turbulence and/or acoustical) 
and surface disturbances. The process by which these 
disturbances are internalized to generate boundary- 
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layer instability waves is referred to as receptivity . 1 A 
fundamental understanding of receptivity is essential 
since it provides the initial conditions for the subse- 
quent growth of the unstable modes that may ulti- 
mately lead to transition. Because of the significant 
wavelength mismatch between the long-wavelength ex- 
ternal disturbances and the much shorter wavelength 
Tollmien-Schlichting (T-S) waves, a length-scale con- 
version must take place to facilitate the energy trans- 
fer. It has been shown theoretically by Goldstein et 
a/., 2- ' 1 Zavobskii et al.J * and Ruban 6 that wavelength- 
reduction mechanisms can range from direct scattering 
associated with abrupt changes in surface conditions 
to the gradual wavelength reduction due to viscous 

I ) 0 u nd ary- 1 ayer grow t h . 

Two general classes of receptivity regions were rec- 
ognized by Goldstein, I) the leading edge regions where 
the boundary layer is thin and growing rapidly and 

II) regions further downstream where the boundary 
layer is forced to make a rapid adjustment. A common 
feature between both receptivity regions is the impor- 
tance of nonparallel effects and is manifested via the 
disturbance motion being governed by the unsteady 
boundary-layer equations. Class II can be further 
subdivided into localized and non- localized receptivity. 
Localized receptivity results from an interaction be- 
tween the freestream disturbances and steady localized 
dist urbance generated by a surface inhomogeneity (e.g. 
humps, gaps, or suet ion /blowing slots). Non- localized 
receptivity stems from an interaction in boundary- 
layer flow with unsteady freestream disturbances and 
the steady disturbance created over surfaces with ex- 
tended regions of short-scale variations (e.g. waviness, 
distributed roughness, uneven suction, etc:.). 

Choudhari 7 investigated the extent of synchroniza- 
tion between the surface waviness and neutral insta- 
bility motion, where for perfect resonance at Branch 
I (lower branch) a u? = at s j . Here a w and a fs de- 
note the wall wavenumber and instability wavenumber, 
respectively and the subscript I refers to conditions 
at Branch I. He shows that the extent of the reso- 
nance region corresponds to a spatial extent of Ax = 

q / o 

0{xiRj' ) where x denotes streamwise distance and 
R Reynolds number. This implies a synchronized sur- 
face waviness will lead to instability amplitudes which 
are larger than those produced by an isolated surface 
roughness by a factor of 0{R' i ), Choudhari 's simple 
analytical result helped to explain some earlier results 
computed by Crouch 8,9 where he found significantly 
larger values for non- localized receptivity versus local- 
ized receptivity. Crouch & Bertolotti 10 have shown 
non- localized receptivity amplitudes to be large over 
a wide range of surface wavenumber. They also show 
an increase of receptivity with surface roughness three 
dimensionality. 

Numerous acoustic receptivity experiments have 
been performed with 2-D surface roughness (e.g. Aizin 


& Polyakov, 15 Saric et al , 12-1 4 Zhou et a/., 15 Koso- 
rygin et a/., 16 Breuer et a/., 17 Kobayashi et a/., 18,19 
and Wiegel k Wlezien 20 ). All of these studies were 
conducted for localized 2-D receptivity except for the 
study done by Wiegel k Wlezien 20 where they con- 
sidered non-localized receptivity. The experiment, of 
Zhou et al. ir> also considered receptivity of an oblique 
surface roughness. That study was conducted for lo- 
calized oblique receptivity where a single strip of tape 
was inclined at angles up to ip w = 45°. They in- 
ferred from their measurements at a fixed streamwise 
location that receptivity was reduced with increased 
oblique angle. They did not consider the reduced lin- 
ear growth rates associated with increased obliqueness. 

The objectives of this research were to further 
understand the acoustic receptivity process for flow 
over wavy surfaces and to correlate theoretical results 
with experimental results. The approach taken was 
to experimentally investigate acoustic receptivity in 
the presence of surface waviness and the subsequent 
boundary-layer disturbance evolution. Previous ex- 
perimental acoustic receptivity data 20 in the presence 
of 2-D surface waviness demonstrated excellent agree- 
ment with theoretical' 9 receptivity coefficients. In the 
presence of oblique surface waviness, previous exper- 
imental data do not exist and inference made from 
previous experimental results 15 for a localized oblique 
roughness strip contradict theoretical findings. The 
current experimental study with 2-D waviness was 
conducted to (i) validate the current approach with 
respect to previous experimental and theoretical re- 
sults and (ii) extend the measurements into regions of 
nonlinear T-S wave development. The measurements 
were then extended to include, for the first time, exper- 
imental data obtained with oblique surface waviness 
for comparison with theory. 10 Receptivity /stability 
experiments are very sensitive to the state of the mean 
flow and environmental disturbances. As a result, 
extreme care was taken to document the flow and en- 
vironmental conditions to avoid ambiguous results as 
discussed by Nishioka k Morkovin 21 and Saric. 22 

Experimental Details 

Facility and Model 

The experiment was conducted in the 2 Foot by 3 
Foot Low- Speed Wind Tunnel located at NASA Lan- 
gley Research Center. The tunnel is a closed-loop 
type with a 10 : 1 contraction ratio. The test sec- 
tion is 0.91m wide by 0.01m high by 0.1m long. The 
turbulence-reduction devices upstream of the contrac- 
tion consist of a honeycomb followed by four stainless- 
steed screens. Speeds of approximately 45 m/s are at- 
tainable in the test section with measured turbulence 
intensities, u/U.^, of approximately 0.1% in the range 
of 0.1 < / < 400 Hz. The test-section floor and ceiling 
are adjustable to achieve a desired streamwise pres- 
sure gradient. Two motorized traverse stages, one with 
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Fig. 1 Top and side views of flat plate model. The 
shaded region represents the location of the insert. 


streamwise travel of 2.1m and the other with vertical 
travel of 150mm, are located just above the test sec- 
tion ceiling. The quoted accuracies of the streamwise 
(j‘) and vertical (y) traverse stages are ±lGG//m/m 
and ±30//m/m, respectively. A 3.8-c??? streamwise slot 
(covered with rectangular wool felt strips to minimize 
air inflow/outflow) along the centerline of the tunnel 
ceiling is provided to accommodate a probe support 
for 2-D traverse motion. 

The model tested was a 12. 7-mm- thick flat alu- 
minum jig plate with a 24 : 1 elliptical leading edge. 
A sketch of the model is shown in Figure 1. Because 
of the potential sensitivity of the receptivity /stability 
characteristics to surface roughness, the upper surface 
(measurement; side) of the model was polished to a 
0.2-jj.m root- mean-square (rms) surface finish. The 
boundary layer on the lower surface of the model was 
tripped using sandpaper grit to establish turbulent 
flow on that surface. The plate w T as equipped with a 
rectangular hole located 24. Gem from the model lead- 
ing edge for the installation/ removal of plate inserts. 
Five inserts were available for the test. One was used 
as the baseline (smooth sample) and four were used 
to mount various samples of receptivity sites. Further 
details of the experimental setup are given by King. 23 

Surface Roughness and Acoustic Excitation 

Receptivity was facilitated by the use of surface 
waviness of height Ah and freestream acoustic exci- 
tation with amplitude u ac . In previous experiments, 
receptivity sites were usually created by applying tape 
or some other adhesive material to the surface of the 
model. For the case of non-Iocalized receptivity, this 
is a very tedious process (e.g. slight misalignment 
angular and/or spatial — of the tape, Ah nonunifor- 
niitios, etc.). An alternate approach for generating 
the surface roughness was considered, which involved 
the use of copper-plated circuit boards. The rough- 
ness patterns were generated using a photolithographic 
process where a high degree of spatial accuracy was 
maintained. In this process a pattern w~as transferred 


Nvh 

Fig. 2 Top and side views of receptivity sites in- 
vestigated — tyu, = 0 for 2-D and i/ j u , ^ 0 for oblique 
roughness. Flow is from top to bottom. 


Sample # 

A h(ftm) 

KuA mm ) 

0 U . (</(<?) 

1 

17.8 

50.25 

0 

2 

35. G 

50.25 

0 

3 

71.1 

50.25 

0 

4 

35. G 

58.88 

0 

5 

35. G 

53.34 

0 

G 

35. G 

50.80 

0 

7 

35. G 

49.53 

0 

8 

35. G 

48. 2G 

0 

9 

35. G 

45.72 

0 

10 

35. G 

50.75 

15 

11 

35. G 

52. 3G 

30 

12 

35.¥ 

55.48 

45 

13 

35. G 

5G.G7 

32 

14 

35. G “ 

54.50 

31 

15 

35. G 

50.27 

29 

10 

35. G 

48.22 

28 


Table 1 Table of roughness samples tested. 



from a photomask to a material layer to produce a 
two-layer pattern with the roughness height Ah be- 
ing equal to the copper plating thickness. The two 
types of patterns considered are represented in Fig- 
ure 2 (the shaded areas represent copper and flow is 
from top to bottom). The first was a 2-D wavinoss pat- 
tern (ip w = 0) with the wavenumber vector parallel to 
the stream direction. The second was an oblique wavi- 
ness pattern where the wavenumber vector makes an 
angle with the stream direction. Highly localized 
disturbances are possible for the oblique patterns as a 
result of the sharp steps formed by the copper strips at 
the leading and trailing edges of the samples. To miti- 
gate this effect, the copper strips were feathered down 
to the substrate material G.4mm from both the lead- 
ing and trailing edges of the samples. The roughness 
patterns considered are presented in Table 1. 

The tunnel was equipped with an array of five 
203.3-mm-diameter woofers (four flush mounted on 
the upstream wall and the other on the downstream 
wall) for acoustic excitation. The speakers were lo- 
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cated as described to minimize unwanted flow distur- 
bances in the tunnel. The source of the controlled 
acoustic field was generated by using a dual-channel 
universal source and amplified with a stereo power am- 
plifier. One channel of the source was used to drive the 
four upstream speakers and the second channel for the 
one downstream speaker. The ability to control the 
absolute amplitude and relative phase of the two chan- 
nels made it possible to generate the required acoustic 
field in the test section. 

Data Reduction 

A single hot-wire probe was utilized to measure 
the streamwise velocity component that represents 
the bulk of the experimental data. A constant- 
temperature anemometer was configured to use a 1 : 20 
bridge with a wire resistive overheat ratio of 1.8. A 
fourth-order polynomial curve fit was found to work 
well for the hot-wire calibration, velocity versus output 
voltage. Broadband and phase-locked velocity fluctu- 
ations were acquired using a typical sampling rate of 
1 kHz. The mean tunnel properties (pressures, tem- 
perature and relative humidity) were measured and 
recorded at each data point. Estimated uncertainties 
of measured mean velocity, E/, are in the range ±1.4% 
for U rU ax ~ 11.1 m/s to ±9.0% for U m i n ^ 1 m/s. Un- 
certainties of nils fluctuating velocity, u y and Reynolds 
number, Hr. are ±4.7% and ±4.0%, respectively. De- 
tails of the error analysis are given by King. 23 

Because of the type of acoustic excitation (i.e. con- 
tinuous wave), the measured rms narrow-band veloci- 
ties UfC 1< * >t included the Stokes wave, the T-S response, 
and any other extraneous response (e.g. probe vibra- 
tion) all at the same frequency f Q . A Stokes w r ave is 
present whenever an oscillating velocity disturbance 
exists in the freestream. Wlezicn et al 2 1,213 found that 
consistent estimates of the T-S modes are achieved 
when the Stokes wave was directly estimated from 
the measured profiles. The various decomposition ap- 
proaches exploit the fact that a large mismatch in 
length scales exists between the Stokes andT-S waves, 
a 95 : 1 ratio for this study. One such approach is 
to make short streamwise surveys that cover approx- 
imately one T-S wavelength, A* s , at a fixed height in 
the boundary layer. 20,26 Over one A ts . the phase of 
the Stokes wave (and for that matter any probe vi- 
bration) is for all practical purposes constant. The 
measured response can then be plotted on the polar 
complex plane. Since the T-S amplitude does not vary 
significantly over one A/*, the centroids (Stokes wave 
± vibration) and the average radii (T-S amplitudes) 
of the off-centered spirals can be computed. The T- 
S wave measured here includes receptivity contribu- 
tions from all possible sources, the leading edge and 
receptivity sites. Such surveys were obtained at wall- 
normal locations of lower maxima of the T-S mode 
shapes. Because the T-S wave is a traveling wave (i.e. 



Fig. 3 Graphical representation of decomposition 
method used for wall-normal surveys. 

linear phase <f> is w.r.t. x), linear regression was em- 
ployed to compute the streamwise T-S wavenumber, 
Of s = \d<pt b /dx\, which was used to obtain wave kine- 
matic results (e.g. \ ts , x = 2tt /a ts . c ph = 2tt f a /a is < 
etc.). 

Most of the disturbance measurements were ac- 
quired by making wall-normal surveys, complete and 
partial boundary-layer profiles, at selected x locations. 
The approach taken here to extract the T-S wave com- 
ponent is a slight variation from earlier methods. 26,26 
Two sets of wall-normal disturbance profiles were ob- 
tained at the same x locations and with the same 
acoustic forcing levels. One set was taken with the 
smooth surface (denoted with ‘'ss" subscripts) and 
the second with a rough surface (denoted by subscript 
“rs”). The disturbance measurements obtained on the 
smooth surface include components of the Stokes wave, 
T-S wave due to leading-edge receptivity, and extra- 
neous disturbances. The disturbances acquired with 
receptivity sites include the same components mea- 
sured on the smooth surface in addition to the T-S 
wave component due to the controlled roughness. This 
is shown graphically in Figure 3. The T-S wave compo- 
nent due to roughness was obtained using the following 

= u Ts e l - u ss e’ <t ’" . ( 1 ) 

Here the T-S component includes receptivity due to 
acoustic scattering at the roughness sites of (1) the 
freestream acoustic field and (2) the leading-edge gen- 
erated T-S wave. The receptivity due to the scattering 
of the leading-edgi* generated T-S wave was assumed 
negligibly small. This was supported by leading-edge 
receptivity measurements in the present study and a 
previous experiment 24 for a 24 : 1 elliptic leading edge, 
both of which indicated small T-S responses. 

Linear stability results presented here were com- 
puted using a computer code by Malik 2 ‘ that incor- 
porates quasi-parallel spatial stability theory. These 
results were used to facilitate the computation of the 
receptivity coefficients, C (= u tii j/u ac ), referenced to 
the Branch I location. Here u ts j denotes the max- 
imum rms T-S amplitude at Branch I and u ac the 
rms freestream acoustic velocity fluctuation. Mea- 
surements of u ts were made downstream of Branch 
I to take advantage of the linear amplification rc- 
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Fig. 4 Neutral stability curve for a Blasius bound- 
ary layer with 2-D disturbances. 

gion for improved signal- to- noise ratios. The values 
of ut#j were obtained by applying linear stability the- 
ory (LST) to the measured values of u ts such that 

C = u ts e- N /u ac (2) 

where N is the amplification factor obtained from 
LST. 

Experimental Approach 

This experiment was conducted with a nominal 
freest ream velocity of U 0 0 = 11.1 m/s and tempera- 
ture of 21 °C. A nominally zero-streamwise-pressure- 
gradient boundary layer was obtained for flow over 
a flat plate (Blasius flow). A planar downstream- 
traveling acoustic wave was used to temporally excite 
the flow with frequency f Q . A simple approach was 
employed to generate the wave train by conducting 
two stream wise surveys in the freest ream above the 
model, each at the same spatial coordinates. One sur- 
vey was taken with only the upstream set of speakers 
activated and the second survey with only the down- 
stream speaker energized. Applying the method of 
superposition, the combined wave field with both sets 
of speakers active was simulated. A scale factor and 
phase shift were selected to give the desired response, 
namely, a constant rms amplitude distribution. 

The' dimensionless frequency of the T-S wave for all 
cases was F — 55 x 10 -6 , which translates to physi- 
cal acoustic forcing frequency of f 0 = 71 Hz. Figure 4 
shows a plot of a neutral stability curve for a Blasius 
boundary layer with 2-D disturbances. A line rep- 
resenting the experimental frequency along with the 
roughness location and extent (At = 0 . 56 / 77 ) is shown 
in the figure. Note that the location of the surface 
waviness was selected such that its midpoint was near 
the location of Branch I to maximize the receptivity. 
The wavenumber of the surface roughness, k w , was 



A/m 

A h- 2 = 2Afn 

> 

w 

II 

it- 

> 

Cl 

fiA/ij 

2e.\Ah i 

4c i A h\ 

C2=2ei 

2e,A hi 

4fi A h\ 

8ei A hi 

e 3 =4fi 

4fi A/ji 

Ah j 

lCfiA/ii 


Table 2 Rotighness heights and forcing levels ex- 
amined — Ahi — 17.8//-T71, ci = 7.6xl0 -5 ( SPL = 84.8). 


selected such that it matched the T-S wavenumber, 
k ts , at Branch I to provide a near- resonance condi- 
tion. Three roughness heights, A //, and three acoustic 
forcing levels, e (= Uac/U^), were examined for the 
2-D waviness. The matrix configuration used to com- 
bine Ah and e are given in Table 2. Detuning was 
examined for the 2-D and oblique roughness patterns. 
This was achieved by maintaining the same forcing fre- 
quency, / 0 , and varying the streamwise wall roughness 
wavenumber, n Uf , while the spanwise wavenumber, 3 W , 
remained unchanged. 

Results and Discussion 

Base Flow Measurements 

Mean Flow 

Normalized streamwise velocity measurements just 
outside the boundary layer, y as 25.4?7?m, revealed 
a pressure recovery region near the leading edge of 
the model with deviations of less than 0.25% of 
from Branch I location to the end of the survey re- 
gion (0.5 < x < 2.03m). The velocities along the 
plate insert location for roughness placement (24. G < 
x < 80.8cm) were all within 1% of Ur^. For the prop- 
erly configured model, the boundary layer developed 
into a nominally zero- pressure-gradient boundary layer 
downstream of the recovery region. For proper com- 
parisons with theory, a virtual origin was defined and 
streamwise locations referenced to this origin are de- 
noted by x v . The virtual origin was computed by first 
acquiring detailed boundary-layer velocity profiles on 
the smooth (baseline) model every 254mm starting at 
x — 508mm. A linear- regression curve fit was applied 
to the square of the measured displacement thickness 
(6*) versus .r. The curve fit was used to extrapolate 
upstream to 5* = 0 and the corresponding x value was 
used as the virtual origin, x = 65. 8mm. 

A more stringent characterization of the flow was to 
the compute the shape factor H for each profile. The 
values of the experimental shape factors were H = 
2.60 ±1.5% in excellent agreement with the theoretical 
Blasius value of H = 2.59. A plot of the measured 
normalize velocity profiles and the theoretical Blasius 
profile versus the normalized wall-normal coordinate i) 
is presented in Figure 5. Excellent agreement between 
the profiles is evident. The two dimensionality of the 
mean flow outside the sidewall contamination region 
(see Figure 1) was validated by obtaining boundary- 
layer measurements at various spanwise locations. 
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u/u„ 


Fig. 5 Normalized velocity profiles at selected 
streamwise locations. 

Freestream Characteristics 

Freest ream turbulent intensity measurements 

{u/U) 0 o, where u is the broadband rms value of 
velocity fluctuations (0.1 </ < 400/fc), were 
obtained using a single hot-wire probe. Values of 
(u/U) oo ~ 0.1% were measured over the entire 

streamwise survey region. More important than the 
magnitude of (u/U)r K , the power spectral density 
of the fluctuating velocity, PSD U . provided relevant 
information on the frequency content of the distur- 
bance environment. Measured PSD U ' s indicate a 
—5/3 power law ( PSD U oc / _5 ^ 3 ), supporting the 
hypothesis that the freest ream disturbances were 
composed predominantly of decaying turbulence. 
Most of the freestrcam energy were contained at 
very low frequencies, below 1 Hz. More than 90% of 
the broadband energy was contained between OAHz 
and 1 Hz illustrating the importance of specifying 
frequency bandwidths when quoting turbulent levels. 
Limited dynamic pressure measurements were ob- 
tained from four pressure transducers flush mounted 
on the test section sidewalls. 23 Correlations of the dy- 
namic pressure signals confirm the absence of coherent 
disturbances in the freestrcam environment near the 
designed forcing frequency, of the experiment that 
can distort the interpretation of the results. 

As a final flow-quality check, the transition Reynolds 
number, Re tr , was measured on the smooth plate. 
This was done by placing the hot-wire probe at a fixed 
location in the boundary layer and slowly sweeping 
through the tunnel speed range. An average value 
of Rctr = 3.2 x 10 r> was measured, indicating low 
freest ream turbulent levels and a smooth model sur- 
face (for comparison, refer to data of Re tr versus 
freest ream turbulence cited in White, 28 pp. 433-436). 
For a configuration with 2-D waviness (Ah = 71. 1/cm) 
and e = 0, a 25% reduction in Rc tr was observed. 
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Fig. 6 Measured freestream acoustic field gener- 
ated in the test section above the model surface — 
(a) amplitude and (b) phase streamwise distribu- 
tions. 

Acoustic- Field Characterization 

Measurements of the controlled freestream acoustic 
field are shown in Figure 6 in the form of (a) am- 
plitude and (b) phase streamwise distributions. Data 
with only the upstream speakers and only the down- 
stream speaker activated are also presented for com- 
pleteness in the figure. In the vicinity of the insert 
region (24.6 < x < 80.8r??i) 5 the variation of the 
acoustic amplitude was small. The phase distribu- 
tion in part (b) of the figure indicates a predominantly 
downstream-traveling wave. Values of the acoustic 
rms amplitudes, u oc , used for normalization in sub- 
sequent sections are the average values measured over 
the plate insert. Streamwise measurements taken at 
various wall- normal and span wise locations verified the 
planar structure of the acoustic field. 
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{deg} 

R, 

At,, /,*{«»»} 

0 

585 

50.2 

15 

583 

50.8 

30 

578 

52.4 

45 

574 

55.5 


Table 3 Tabulated data of Reynolds numbers and 
streamwise T-S wavelength at Branch I for distur- 
bance wave angle 0^ obtained from LST. Note that 
R % 85.92 * x v where x v is in centimeters. 

Receptivity Measurements 

Effect of e, Ah, and 0 U . on Receptivity 

The influences of acoustic forcing level, e , roughness 
height, Ah, and roughness obliqueness, ij > w , were ex- 
amined by conducting measurements with the matrix 
arrangement shown in Table 2 and utilizing samples 
1—3 and 10 — 12 of Table 1 . At the location of Branch 
I (set 1 Table 3), the largest roughness height, Ah 3 , nor- 
malized by 5* was 0.05 and the corresponding rough- 
ness Reynolds number was RcaHj = 47. Therefore, 
all roughness heights considered were well within the 
roughness Reynolds number range of a hydraulically 
smooth surface (see White, 28 pp. 436-439). Measure- 
ments were obtained over a range of Reynolds num- 
bers, R > 570 (x > 50.8cm), that extended beyond 
the location of Branch II. The results to be presented 
here were obtained for the condition where the wave- 
length of the wall roughness X w matched that of the 
T-S wave at Branch I, the near-resonance condition. 
In this context, the term near-resonance describes the 
condition when = c\ ts j which does not account for 
the finite speed of the freestream acoustic wave, i.e. 
o«c / 0 . 

Measurements of T-S amplitudes at locations just 
downstream of Branch I indicated a virtually linear 
response of u ls with both u ac and Ah in the presence 
of 2 -D roughness. T-S amplitude eigenfunctions are 
shown in Figure 7 at R — 1038 for various forcing 
products, f • Ah. The rms T-S amplitude, u ts , w T as 
normalized by u ac and h where h is the step height 
normalized by the smallest step) height. The linear 
behavior is evident from the data collapse in the fig- 
ure and indicates linear receptivity for all acoustic and 
wall forcing levels considered. The 180° phase shift 
at rj ~ 25* is evident in the plots of the T-S phase 
(pts (not shown). Measured data depicting the stream- 
wise evolution of the T-S eigenfunctions indicate good 
agreement with LST for locations upstream of where 
nonlinear effects first became apparent in the near- 
wall region. The nonlinear effects became evident at 
the same x location for a given value of e- Ah, indepen- 
dent of the individual values of e or Ah. More on the 
nonlinear evolution of the T-S instability is discussed 
later. 

To investigate the influence of T-S wave oblique- 
ness on receptivity, samples 2 and 10—12 (refer to 
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Fig. 7 Normalized T-S amplitude profiles with 2-D 
surface waviness for various values of e and Ah. 



Fig. 8 T-S amplitude profiles for various values of 
0 U . at R — 1038. 

Table 1) were used to generate the oblique T-S waves 
by the application of oblique waviness. The roughness 
wavenumbers were selected to create a near-resonance 
condition at Brancli I for each 0 U , considered (see Ta- 
ble 3 for tabulated values of the streamwise T-S wave- 
length at Branch I, A t s j,x (= 27t /c\i s j)). Measured 
and LST T-S amplitudes at R = 1038 (x = 152.4cm) 
are presented in Figure 8 for various y> w . The LST 
results were normalized to match the experimental am- 
plitudes at the maximum of u ts for these profiles. Note 
that the relative peak (inner and outer) amplitudes of 
the measured T-S profiles do not scale directly with 
LST results. This may be due in part to nonparallel 
effects not accounted for in LST since nonparallelism 
becomes more important for three-dimensional distur- 
bances (see Bertolotti 29 ). Nevertheless, the salient 
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Fig. 9 Streamwise amplification data for various 
values of ij > w . 

features of the eigenfunction profiles were captured by 
LST. For 2-D disturbances, a single peak in the near- 
wall region below the critical layer ( 7 / ^ 1) was evident. 
As the dist urbances became oblique, some residual of 
this peak remained and the appearance of a dual peak 
structure in the near-wall region — which straddles 
the critical layer — is evident. An upward shift of 
the global maximum to a wall-normal location above 
the critical layer was observed. The amplitude of the 
apparent residual peak reduced with increased oblique- 
ness. An upward shift of the wall-normal location of 
the 180° phase jump was also evident. 

The streamwise amplitude distribution of the mea- 
sured T-S eigenfunctions were compared to the am- 
plification predicted by LST. Figure 9 shows a plot 
of JV-factor from LST versus Reynolds number along 
with the experimentally measured maximum T-S am- 
plitudes. The experimental amplitudes (ln(u ts /u ac )) 
were obtained from measured eigenfunction profiles 
taken at selec ted .7' locations and scaled for comparison 
with LST. Excellent agreement between the measured 
and LST data for 2-D waviness is evident except for 
the large values of u ts measured for /? ~ 717 just- 
downstream of Branch I (in a region over the wavy 
surface). The' large values of u tif measured in the near- 
wall region just downstream of Branch I may partly be 
explained by the non-monotonic increase of the total 
instability response resulting from the localized peak 
(overshoot) in the receptivity region clue to the exter- 
nal forcing. ‘T 9 The location of Branch II is in excellent 
agreement- with LST. The data also demonstrate (not 
shown) that consistent amplification curves are obtain- 
able by tracking either the inner or outer peaks, pro- 
vided the response is linear. Excellent agreement was 
also observed between the measured and LST amplifi- 
cation curves for small oblique angles. The agreement 
became less attractive for the larger oblique angles, 



Fig. 10 Non-locaiized 2-D receptivity coefficients 
from present experiment (F = 55 -X 10 r> ) and past 
experiment (F % 46 x I0~ 6 ) as a function of Ah. 

particularly near the location of Branch II. Again, this 
may be a result of nonparallel effects not captured by 
LST. Agreement between the measured and experi- 
mental amplifications is important since theory (e.g. 
LST) is required to backtrack the effective Branch I 
amplitudes in an accurate fashion. 

The receptivity coefficients, C, to be presented were 
referenced to Branch I and computed utilizing Equa- 
tion 2. The stipulation that all values of Ut s employed 
were obtained downstream of Branch I and the rough- 
ness surface and upstream of Branch II was imple- 
mented. Figure 10 shows a plot of the experimental 
2 -D receptivity coefficients (plus linear regression fit) 
versus the normalized step height. Data from the ex- 
periment of Wiegel Wlezien 20 for a different nondi- 
mensional frequency (F & 46 x 10 -6 ) are also included. 
The symbols represent mean values and the error bars 
the maximum plus/minus deviations observed in the 
measurements for all values of c. The leading-edge re- 
ceptivity coefficient for the smooth surface (/i = 0) 
was calculated using data obtained in short stream- 
wise surveys. The small values of the leading-edge 
receptivity coefficients (C « 0,035) support the de- 
composition assumption discussed previously. Good 
agreement is observed between the present and previ- 
ous 20 experimental receptivity coefficients. The cor- 
responding T-S wave kinematics were inferred from 
data obtained with short streamwise surveys (i.e. 4>t s 
versus x data) where (= -depts/dx) was com- 
puted directly. The average T-S wave phase speed, 
Cph/Uoo = 0.33 ( c p h = 2Tzf 0 fa ts ), at R = 1038 for all 
values of e and Ah considered is in excellent agreement 
with LST result. c. v \ x j = 0.33. 

Figure 11 shows a plot of measured receptivity co- 
efficients as a function of rp w along with theoretical 
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Fig. 11 Measured (F — 55 x 10 6 ) and theoretical 
(F — 56 x 10 _f5 ) receptivity coefficients versus ij.i w • 

data (F = 5G x 10~°) by Crouch k Bertolotti. 10 The 
agreement between the measured and theoretical data 
is good, thereby providing experimental evidence of 
increased receptivity with increased wave obliqueness. 
The confidence level for the value of the experimental 
receptivity coefficient at i}) w = 45° is not very high 
since only one data point was used to obtain the value 
of C. This may be one explanation for the large dif- 
ference observed there along with the fact that the 
mode- interaction theory used to compute the theoret- 
ical values of C degrades with three-dimensionality, 
which is largely attributable to nonparallel effects. 10 
The overall agreement between the measured and LST 
Cph/Utx, with respect to ij) w is good, within the exper- 
imental uncertainty of the measurements. 

Effect of detuning on Receptivity 
In this section, results are presented for the so- 
called detuned conditions. Detuning is referred to as 
the condition when there is a mismatch between the 
wavelengths of the wall roughness and the T-S mode 
(\ tv ^ \( S ) for a given disturbance frequency. The en- 
ergy transfer from the forced disturbance becomes less 
efficient as the degree of detuning is increased. As a 
measure of the detuning, the detuning parameter a is 
defined as 


In previous studies (e.g. Wiegel k Wlozien 20 ), detun- 
ing was generally examined by varying C/ 00> which 
consequently varied a*.*,/. This was done because of 
difficulties realized in varying a w . Because of the type 
of receptivity sites used in the current study, a was 
varied by changing <y w . Detuning data for both 2- 
D and oblique waviness were obtained using samples 
2 and 4-0 and samples 11 and 13 — 1G, respectively 
(refer to Table 1). All the data were obtained with the 


Fig. 12 Measured (F = 55 x 10~ n ) and theoreti- 
cal (F = 55 x 10 -fi and F — 5G x 10 -6 ) receptivity 
coefficients versus the detuning parameter. 

forcing product 62 ■ A//j. 

The maximum receptivity was not expected at the 
near-resonance condition a llf = a*#,/ as predicted by 
receptivity theory that assumes incompressible flow 
(infinite acoustic wave speed). To account for a finite 
acoustic wave speed, one lias to take into consideration 
the wavenumber of the freestream acoustic wave, o rtC , 
as conjectured by Wiegel k Wlezien. 20 For the present 
experiment, care was taken to acoustically excite the 
freestream with a downstream- traveling wave. Clioud- 
hari k Street! 30 discussed the finite wave speed effect 
and showed that a resonance condition exists when 
(c*w,M = {atsj - <XaciPt8,i) since ;3 ac = 0 here. 

Measured and theoretical 7,0, 10 receptivity coeffi- 
cients are presented in Figure 12 for various values 
of the detuning parameter (modified) a'. The rough- 
ness height used to compute the theoretical receptivity 
coefficients 7,9, 10 for single-mode waviness was the first 
Fourier coefficient of a square wave pattern with height 
A/i, namely the mode corresponding to the funda- 
mental wavenumber, a w < with amplitude Ah/n. The 
parameter a f shifts the values of a in Equation 3 to ac- 
count for the finite acoustic wave speed. The detuning 
parameters are related by a f — a + ci ac /&tsj where the 
respective values of o GC /a ts j for the 2-D and oblique 
waviness are 1.04% and 1.08%.. The receptivity coeffi- 
cients are normalized by C a where C a is the maximum 
value of C (measured or computed) except for the the- 
oretical C a value for the oblique roughness which is the 
theoretical value of C at the same a' (~ —2.8) as the 
corresponding experimental value. Good agreement is 
obtained between the measured and theoretical values 
demonstrating the effectiveness of non- localized recep- 
tivity over a wide range of a w . A maximum difference 
of less than 10% was obtained between the measured 
maximum and the corresponding theoretical value of 
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Fig. 13 Streamwise amplification data with 2-D 
disturbances for various values of t ■ Ah. 

the receptivity coefficients, C a . 

Disturbance Evolution 

In this section, data describing the evolution of the 
boundary-layer disturbance field are presented. Em- 
phasis is placed on results obtained in the presence of 
2-D surface roughness since nonlinear behavior was ev- 
ident only for large values of the forcing combination 
(f • Ah). The broadband results obtained with oblique 
waviness were qualitatively similar to those obtained 
with 2-D waviness and intermediate forcing, C2 * A/cj. 
For this reason, oblique results are not discussed ex- 
plicitly in this section except where minor differences 
were noted. The measured streamwise amplitude dis- 
tribution of the normalized T-S eigenfunctions were 
compared for various values of e* Ah in Figure 13. The 
solid symbols represent data of the outer peak and the 
open symbols data of the inner (maximum) peak. Non- 
linear effects are evident from the amplification curves 
for the forcing products of * A/&3 and 63 * A/e* (see for 
comparison Figure 9). For the case of €2 * A/13, non- 
linear effects are apparent only in the near-wall region 
and first become discernible when u ts /U 0 0 > 0.8% at 
7? ~ 1120 (just upstream of Branch II, Rjj — 1129). 
For the largest forcing product (c 3 • AA3), nonlinear- 
ity is evident in both the near-wall and outer regions 
of the boundary layer. The nonlinearity was observed 
first near the surface and eventually spread through- 
out the shear layer due to the interaction of normal 
vortidty, which is driven by the mean shear (see Co- 
hen, Brener, k Haritonidis 31 ). Nonlinearity was first 
observed for Uts/Uoc > 1.0% at R ~ 1034 ( R < 7?//). 

Broadband (2 < / < 4i)0Hz) rms fluctuating veloci- 
ties, u, measured for e 2 A h 2 are presented in Figure 14. 
For this intermediate forcing level, the flow remained 
laminar (see Figure 13). In the early stages (R < 943, 
x < 127 cm), the rms velocity profiles resemble those 

10 



Fig. 14 Measured broadband rms fluctuating ve- 
locity for forcing combination of <2 * A h >. 

of Klebanoff modes as measured by Klebanoff, 32 and 
Kendall. 33,34 The profiles are indistinguishable from 
profiles measured on the smooth model surface with 
no controlled acoustic forcing (baseline configuration) 
— e — 0, Ah - 0. The maximum rms fluctuations 
occur near the midpoint of the boundary-layer thick- 
ness ( 7 } ~ 2.3). The rate of increase of the maximum 
rms velocity at each streamwise station is small for 
these low- freest ream turbulent intensity levels. For 
measurements taken further downstream (7? > 943), 
the contributions from the T-S wave become dominant 
in the near-wall region and near the boundary- layer 
edge due to the dual-lobe shape of the T-S eigenfunc- 
tions. Near the midpoint of the boundary layer, the 
rms velocities are still consistent with those measured 
on the baseline configuration. As shown using parab- 
olized stability theory by Herbert and Lin, 35 the T-S 
waves and these low-frequency Klebanoff modes ap- 
pear to exist independently, but both may contribute 
to the intensity at a given frequency. For 63 • A/13 
where nonlinear effects and consequent breakdown oc- 
curred, increases of one to two orders of magnitude in 
the maximum values of u were observed. Power spec- 
tral densities of the fluctuating streamwise velocity, 
PSD U , corresponding to data presented in Figure 14 
are shown in Figure 15. The spectra are shown for data 
obtained at ?/ ~ 2.3. The controlled acoustic forcing is 
clearly evident in the spectra at f 0 = 71 Hz. The spec- 
tra shown consist of significant low-frequency energy, 
a characteristic of Klebanoff modes. Modest amplifi- 
cation at frequencies (/ < 2077T) below the unstable 
T-S modes is evident at R = 1205 (x = 203.2cm) and 

was also observed for the baseline configuration. As 

expected, significantly more amplification is evident 
in the unstable T-S band. 
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Fig. 15 Spectra of fluctuating velocity near r] % 2.3 
for forcing product of f 2 A/* 2 . 

To investigate the coherent nature of the boundary- 
layer disturbance field, particular consideration was 
givcai to disturbances at the fundamental (/ 0 ), sub- 
harmonic (f 0 / 2), and first harmonic (2 f Q ) frequencies. 
The coherence function, 7, (between the designated 
instability wave and the input signal used for the 
acoustic forcing) at the three frequencies of interest 
was computed. Significant values of 7 for 62 * A/12 were 
obtained at the fundamental and harmonic frequen- 
cies for a wall-normal location corresponding to the 
maximum of the T-S wave amplitude; however, negligi- 
ble coherence was observed for the subharmonic band. 
The amplification curves at the selected frequencies 
reveal that the harmonic disturbance essentially mir- 
rors the amplification of the fundamental disturbance. 

As a result, the harmonic disturbance was believed 
to be merely driven by the fundamental and a result 
of a nonlinear self- interaction. No amplification was 
observed at the subharmonic frequency. Coherence 
and amplification data for C3 ■ A//3 are shown in Fig- 
ure 16. A sharp roll off in coherence of the fundamental 
and harmonic disturbances are evident for R > 1038 
(x > 152.4cm) indicating the onset of significantly 
nonlinear behavior leading to turbulence. This is in 
agreement with results inferred from the amplification 
curves in Figure 13. Again, negligible coherence was 
observed at the subharmonic frequency. Referring to 
part (b) of the figure, it is evident from the amplifica- 
tion curves that the harmonic disturbance grew at rate 
different from that of fundamental, indicating a strong 
nonlinear interaction at the harmonic frequency. The 
growth seen at the sub harmonic frequency was not due 
to phase-locked amplification, but a result of broaden- 
ing of the spectra indicative of the breakdown process. 
The transition Reynolds number for this scenario was 
Re tr ^ 1.2 x 10 6 which is approximately one-third of 
the baseline Re tr (~ 3.2 x 10 6 ). 
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Fig. 16 Data showing the fundamental, subhar- 
monic, and first harmonic (a) coherence functions 
and (b) growth curves for €3 • Ah 3. 

Finally, broadband spectra for the intermediate and 
largest forcing levels are presented in Figure 17 for a 
range of stream wise locations (717 < R < 1205 or 
76.2 < x < 203.2r??7) at wall- normal locations corre- 
sponding to the maximum of u ts . For the intermediate 
forcing level (^2 * A/^) in part (a) of the figure, the 
energy and associated growth at the fundamental forc- 
ing frequency (f Q — 71 Hz) and its first harmonic are 
clearly evident. The energy at / = 120 Hz is associ- 
ated with electronic noise as evident from the constant 
streamwise amplit ude. Amplification of broadband en- 
ergy centered about / « GO Hz presumably a result 
of the natural tunnel disturbance environment is dis- 
cernible in the spectra for > 127.0cm {R > 943). 
Broadband amplification is also evident about / « 
130 Hz for x > 152.4cm (R > 1038). This latter 
broadband energy is possibly an interaction between 
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Fig. 17 Spectra of streamwise velocity taken at over a range of streamwise locations for (a) 

C 2 *A /*2 and (b) t^-Ahs. The subharmonic, fundamental, and first harmonic frequencies (along with 120-/7,3 
electronic noise) are denoted in part (a) of the figure. 


the forced T-S wave and the naturally occurring broad- 
band T-S field. The broadband amplification about 
/ ^ 130Hz was not evident for the oblique waviness 
configuration (0V = 30°). This is believed to be due to 
the reduced energy amplitudes of both the forced T-S 
wave and the naturally occurring broadband T-S field. 
The situation was quite different for the largest forcing 
combination in part (b) of the figure. The presence 
of signific ant energy at the fundamental and higher 
harmonic frequencies is clearly discernible at the on- 
set of nonlinearity (j* > 152.4cm or R > 1038). The 
harmonic bands broaden until a developed turbulent 
spectrum was obtained. Broadband amplification in 
the unstable T-S band of frequencies as seen for the 
intermediate forcing level above was not discernible 
in the present spectra. Similar harmonic- dominated 
spectra were obtained in the experiment of Brener et 
ai 17 This type of breakdown From laminar to tur- 
bulent flow is indicative of a K-type breakdown as 
obtained by the classical results of Klebanoff et a /. 30,37 
The selectivity of the single-frequency wave, imposed 
by the controlled acoustic forcing and receptivity sites, 
created a dominant primary T-S wave with amplitudes 
much larger than other two- and three-dimensional 
modes resulting from background disturbances. The 
primary mode at the fundamental frequency domi- 
nated the breakdown process and resulting nonlinear 
interactions produced harmonics as is evident in the 
disturbance spectral evolution (sec Figure 17(b)). 

Concluding Remarks 

The most widely used transition prediction tool is 
currently the e N method, which is a semiempirical 
method based on linear amplification (not disturbance 
amplitudes) of unstable modes. This approach does 
not include any physics regarding the role of the ex- 
ternal disturbance environment on the boundary-layer 


instability development. Improved transition predic- 
tion tools based on an amplitude approach need to 
incorporate physics relating the effects of the environ- 
ment on the transition process. This approach needs, 
as an input, information about the initial boundary- 
layer disturbance amplitude, and this need is met by a 
better understanding of the receptivity process. Con- 
sequently, receptivity plays a crucial role in any tran- 
sition prediction tool that utilizes an amplitude based 
approach. The objective of the current research was 
to better understand the acoustic receptivity process 
and subsequent disturbance evolution in the presence 
of wavy surfaces and correlate results with theory. Ex- 
perimental investigations were conducted to examine 
acoustic receptivity and subsequent boundary-layer in- 
stability evolution in the presence of 2 -D and oblique 
surface waviness. These investigations provided a data 
base which can ultimately be used to validate cur- 
rent and fut ure theoretical receptivity /transit ion mod- 
els needed to improve the current understanding of 
1 am i n ar- to - 1 u r b u lent transit ion . 

A subtle accomplishment of the current study was 
the ability to introduce a variety of steady boundary- 
layer disturbances via the well-defined receptivity sites 
without any measurable first-order effects on the base 
flow. This study explored, in addition to 2-D acoustic 
receptivity, the subsequent evolution of the boundary- 
layer instability. The current acoustic receptivity mea- 
surements with 2-D waviness are in excellent agree- 
ment with the existing experimental measurements 20 
and theoretical results . 7,9 For a configuration where 
laminar-to-turbulent breakdown occurred (^3 * A/ej), 
the breakdown process was found to be dominated by 
energy at the fundamental and higher-harmonic fre- 
quencies, indicative of K-type breakdown. Receptivity 
data were acquired experimentally for the first time 
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over oblique wavy surfaces (0 U , = 15°, 30°, and 45°) 
to measure the dependence of receptivity on roughness 
obliqueness (or three dimensionality). The current 
findings were in substantial agreement with theory 
and experimentally confirmed the earlier theoretical 
results 10 that receptivity increases markedly with wave 
obliqueness. Boundary-layer transition was not ob- 
served for the forcing combination (ez * A// 2 ) used with 
the oblique surface waviness. 

The transition Reynolds number measured on the 
baseline configuration was Rc tr « 3.2 x 10°. The value 
of Rc t j for the 2 -D surface roughness with a forcing 
combination of 63 • A /13 was significantly less (~ 60% 
reduction) than that for the baseline configuration. 
The drastic shifts in the value of Rct r demonstrate the 
sensitivity of the transition location to the freestream 
environment and surface conditions. 
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